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Abstract. We use expansion-normalised variables to investigate the Bianchi 
type VIIo model with a tilted 7-law perfect fluid. We emphasize the late-time 
asymptotic dynamical behaviour of the models and determine their asymptotic 
states. Unlike the other Bianchi models of solvable type, the type VIIo state 
space is unbounded. Consequently we show that, for a general non-inflationary 
perfect fluid, one of the curvature variables diverges at late times, which im- 
plies that the type VIIo model is not asymptotically self-similar to the future. 
Regarding the tilt velocity, we show that for fluids with 7 < 4/3 (which in- 
cludes the important case of dust, 7 = 1) the tilt velocity tends to zero at 
late times, while for a radiation fluid, 7 = 4/3, the fluid is tilted and its vor- 
ticity is dynamically significant at late times. For fluids stiffer than radiation 
(7 > 4/3), the future asymptotic state is an extremely tilted spacetime with 
vorticity. 



1. Introduction 

Much progess has been made recently in understanding spatially homogeneous 
(SH) cosmologies containing a perfect fluid with an equation of state p = (7 — l)p, 
where 7 is a constant. For the SH cosmologies known as the Bianchi models, the 
universe is foliated into space-like hypersurfaces (defined by the group orbits of 
the respective model) and the Einstein equations describe the evolution of these 
hypersurfaces [21 EH El 00 EU ■ For these models there are two naturally defined 
time-like vectors: the unit vector field, n^, normal to the group orbits, and the 
four- velocity, it p , of the perfect fluid. When analysing the Bianchi models it is 
common to utilize as the preferred timelike vector field in the formalism, while 
the fluid velocity may or may not be aligned with If is not aligned with 
71**, the model is called tilted, and non-tilted (or orthogonal) otherwise [7j- 

All of the non-tilted Bianchi models, except for the type IX model 1 , have been 
studied in detail (see e.g., EH)- For tilted Bianchi models the picture is less 
complete; the type II model was studied in JT], the type V model in [T^l ITTA fT~H 
HTn ITT)] , type VI in QH GS|, types IV and VIL, in [BO ED], a subset of the type 
Vl/i models in and irrotational type VIIo models were studied in JHj- In this 
paper we will use a dynamical systems approach and investigate the tilted type 
VIIo models. 

In the class of spatially homogeneous Bianchi models the type VIIo model is a 
special, and particularly interesting, case. For Bianchi models containing a non- 
tilted perfect fluid, it was shown that the type VIIo model experiences a self- 
similarity breaking at late times PH][22] (see also [221 IS])- The irrotational models 
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Figure 1. The Bianchi generalisations of the closed (k = 1), flat 
(fc = 0) and open (k = -1) FRW model. 



show the same behaviour |19j . The reason for this self-similarity breaking is because 
one of the curvature variables grows unbounded, which leads to an oscillation of the 
shear and the curvature variables with a frequency increasing with time. We will 
see that the same behaviour occurs for the fully tilted type VIIo models; however, 
in this case there are two oscillations with different frequencies (one in the shear 
and curvature variables, and the other in the tilt velocity). 

The investigation of the type VIIo model is important for several reasons. First, 
the type VIIo model is the most general spatially homogeneous model allowing for 
the flat Friedmann-Robertson- Walker (FRW) model as a special case (see Fig0. 
Both the Bianchi type I model and the VIIo generalise the flat FRW model; however, 
since the Bianchi type I perfect fluid model does not allow for tilt, only the Bianchi 
type VIIo model can be used to investigate the effect of tilt on the evolution of 
the universe close to flatness. In addition, the type VIIo model plays an important 
role in the Bianchi hierarchy (depicted in Fig|2J). The higher up in the hierarchy, 
the more free parameters the model contains and thus the more 'general' they are. 
Fig. |3also shows the possible limits (solid arrows) of the various Lie algebras. For 
example, the type VIIo is a special limit of both the 'most general' type VIII and 
type IX models. In order to understand a particular Bianchi model, a complete 
understanding of all its descendants at the lower levels is needed. This implies that 
an understanding of the type VIIo model is necessary for a full understanding of 
the two general semi-simple models (i.e., the type VIII and IX models). 



2. Equations of motion for the tilted VIIo model 
The Bianchi type VIIo model has the Lie algebra of Killing vectors given by: 

This Lie algebra is solvable and in ^5] a formalism for studying the tilted solvable 
Bianchi models was given. In this paper we will use this formalism to study the 
Bianchi type VIIo model with full tilt. 

In the dynamical systems approach we introduce expansion-normalised shear 
variables (bold variables are complex variables) (E + ,S x ,Si) and the curvature 
variables (A, N, N x ). More specifically, N a a = 2\/3N while N x consists of the 
trace-free part of N a b- Moreover, for the type VIIo model, A = 0. Regarding the 
fluid, fl is the expansion-normalised energy density and the tilt velocity is given by 
one real and one complex variable, v\ and v = V2 + W3. We note that the fluid in 
the Class A models has non-zero vorticity if and only if N a bV a ^0 0, which implies 
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Figure 2. The Bianchi Hierarchy: Solid arrows indicate Lie al- 
gebra contractions, or certain special limits of the Lie algebra pa- 
rameters. Dashed lines indicate one-parameter families of non- 
isomorphic Lie algebras. 



that the type VIIo model has vorticity 2 if and only if v 7^ 0. It is also necessary to 
introduce the dimensionless time variable, r, defined by 

(1) *=I 
W dr H' 

where t is the cosmological time and H is the Hubble scalar. 

The equations of motion in |19j are given in terms of an arbitrary gauge function 
</>. In terms of the above variables, the gauge transformation is given by 

(2) (N x ,E x ,Ei,v) - (e 2 ^N x ,e 2 ^£ x , e ^£x,e^v). 

Several different gauges were proposed and discussed (including their advantages 
and disadvantages). In this paper we will adopt the F- gauge' for which <p' = 0. 
This still leaves us with an unspecified constant gauge freedom. The system of 
equations will therefore contain one more variable than the number of physical 
degrees of freedom. The physical properties of the system can be extracted by 
considering gauge independent quantites such as, for example, S X N* and N* v 2 . 



2 From 1111 (see Appendix C), the fluid vorticity W a of the Bianchi VIIo model (components 
with respect to the G3-adapted frame) is given by 



Wi 



1 1 



-N CD v c v D v 1 , W A 



1 

2B 



2B1-V 2 

where the upper case subscripts run from 2 to 3, and 



Na°vc + - — \j^N cd v c v d va 
1 — V z 



Wo = -v a W a 



B 



G-Vl - V 2 
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In the 'F-gauge' (</)' = 0) the equations of motion are: 

(q - 2)£ + + 3|S 2 | 2 - 2|N X | 2 + ^- {-2v\ - 

)£ x + - 2V3NX7V + ^^v 2 

2G+ 



(3) 


K 


= (9 


(4) 


K 


= (9 


(5) 


si 


- (<z 


(6) 


N' x 


= (<Z 


(7) 


N' 


= (<Z 



%/3 7 ft 



Vl 



The equations for the fluid are 



(8) O' = — {2<z-(3 7 -2) + [2g( 7 -l)-(2- 7 )- 7< S]F 2 } 

(9) v[ = (T + 2S + )w! - 2%/3Re(Siv*) + V3Im(N* x v 2 ) 

(10) v' = (T-T, + -iV3Nv 1 )v-V3(E x +iN K v 1 )v* 



v(i-v 2 ) 

where 



q = 


2E2+ l(3 7 -2) + (2 7 )v 2 

2 i + ( 7 -i)v 2 


s 2 = 


S 2 h + |S x | 2 + |S 1 | 2 


5 = 


S ab c a c fc , c a c Q = 1, v a = Vc a 


V 2 = 


vl + W\ 2 , 


G+ = 


l + ( 7 -l)^ 2 , 


T = 


(3 7 -4)(l- V 2 ) + (2- 7 )V 2 S 


l-( 7 -l)F 2 



(12) 

These variables are subject to the constraints 

(13) 1 = £ 2 + |N x | 2 + r! 

(14) = 2Im(S* x N ^ ' lQVl 



(15) = iSi^ + iS*N> 



G+ 

+ o7 



The parameter 7 will be assumed to be in the interval 7 € (0,2). The gener- 
alised Friedmann equation, ea. (|13fl . yields an expression which effectively deter- 
mines the energy density ft. The state vector will therefore be considered to be 
X = [£+, S x , Si, N X) N, Vi, v] modulo the constraint equations 1)14(1 and l)15|l. 
Thus the real dimension of the dynamical system is eight of which one is the re- 
maining gauge freedom; i.e., the dimension of the physical state space is seven. 
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The dynamical system is invariant under the following discrete symmetries: 

[E+,£ x ,£i,N x ,JV,t;i,v] >-> [£ +) E x ,£x,-N x ,-jV,-ui,-v] 
[E+,E x ,Ei,N x ,JV,«i,v] ^ p + ,u 2 r x ,-uE;,u 2 N* x ,]V,- Vl ,jv*], u e S 1 



These discrete symmetries imply that without loss of generality we can restrict the 
variables N > and «i > 0. 3 We also note that the free parameter u in <p2 is, in 
fact, the remaining gauge transformation. 

2.1. Invariant Subspaces. In this analysis we will be concerned with the following 
invariant sets: 

(1) T(V7T ): The general tilted type VII model with N 2 - |N X | 2 > 0. 

(2) Ti(VH ): The irrotational type VII models defined by iV 2 - |N X | 2 > 0, 
Ex =v = 0. 

(3) F{VIIq): The set of fixed points of 02 (without loss of generality we can 
set u = 1) defined by Vi = Im(E x ) = Re(Ex) = Im(N x ) = Im(v) = 0. 

(4) B(VIIq): Non-tilted Bianchi type VII models with N 2 - |N X | 2 > 0, v x = 
v = 0. 

(5) T(II): The general tilted type II model given by N 2 - |N X | 2 = 0. 

(6) B(I): Type I: N = N x = V = 0. 

(7) dT(I): "Tilted" vacuum type I: Q = N = N x = 0. 

We note that the closure of the set T(VIIq) is given by: 



(16) T(VII ) = T(VII ) U T{II) U B{I) U 9T(7). 



Morever, in T(VIIq) we have the bounds 

(17) + |S x | 2 + |Sx| 2 + |N x | 2 <1, w? + |v| 2 <l; 

hence, all variables are bounded except for N (which can become arbitrary large). 

2.2. Monotonic functions. There are three monotonic functions which are useful 
for our analysis: 



(1 - F 2 )K 2 -7) 6 
(18) Zi = afi 1 ^, a = 1 Gl J vr , r = - 7 , 

Z[ = [2(l-r)g+(2 + 2r-37)+T5]Z 1 . 
We note that (by the same trick as in |19|) 

2(i-r) g + (2 + 2r-3 7 ) + r5 

1, w 9 , ,9s 7[l8(l-7) + (10-97)l/ 2 l 
> -(14 - 15 7 )(2E 2 + |N X | 2 ) + 1^ U—Lsi. 

Thus Z\ is monotonically increasing in T(VIIq) for 7 < 14/15. 



^There is a sublety regarding this choice since, in general, v\ = is not an invariant subspace. 
The state space can be considered an orbifold with a mirror symmetry at v\ = 0; in particular, 
this means that any equilibrium point in the region v\ > has an analogous equilibrium point in 
the region v\ < 0. 
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(^-|N x | 2 )"/3n a _ (1-^)1(2-,) i 

(1 + nS+) 2(1+n) G + 4 

7' 4 ( / v , s2, n 2mv |2 , ( 1 + »)( 2 -5n) 2 

^2 = 7— v S ( s + + "J +(l-n S x + Si 

1 + nh + I 2 

+ ( i + 2,,Ki + n) n [2(1 _ ntf + (2 _ 5 „ )|vn | z , 

Z2 is monotonically increasing for the following subspaces: 4 Ti(V77o) f° r < 7 < 
2, and T(VJJ ) for < 7 < 6/5. 



(20) 



[iV|v| 2 +Re(N* x v 2 )] : 



(1_ 72)2(2-7)^ ' 

Z' 3 = 3(5 7 -6)Z 3 . 

Thus is monotonically decreasing (7 < 6/5) or increasing (6/5 < 7) in T(VIIq)\ 
Tx{VII Q ). 

3. Qualitative analysis 

The above monotonic functions allow us to obtain some results regarding the 
asymptotic behaviour of Bianchi type VIIo universes: 

Theorem 3.1. For < 7 < 6/5, all tilted Bianchi models (with > 0, V < I) of 
type VIIq are asymptotically non-tilted at late times. 

Proof. For < 7 < 14/15 we can use the monotonic function Zi, which immediately 
gives the desired result. For 2/3 < 7 < 6/5, and using the monotonic function Z3, 
we get 7V| v| 2 + Re(N x v 2 ) — > 0. Moreover, using Z2, we get N — > 00 and hence 
|v| — > 0. This implies that the solution asymptotically approaches the invariant 
subspace Ti(VJJo). The Ti(VII ) analysis J^I can now be applied which shows 
that v-l — > 0. The theorem now follows. □ 

In fact, we believe that the type VIIo models are asymptotically non-tilted for 
< 7 < 4/3; however, this is not covered by this theorem. 

Another important observation is the divergence of the variable iV at late times. 
In fact, we firmly believe that: 

Conjecture 3.2. For a non-inflationary perfect fluid (2/3 < 7 < 2) and any initial 
condition in T(VIIq) with V < 1 and Q > 0, we have that 

lim IN] =00. 

r — >+oo 



There are several results that support this conjecture. First, use of the monotonic 
function Zi immediately establishes this result for T\{VIIq) and for 2/3 < 7 < 6/5 
with general tilt (which includes the important case 7 = 1). Moreover, using Z3 

^We note that all of these functions are also monotonic for the Vlo model, where the importance 
of the value 7 = 6/5 is more clear. In fact, using these monotonic functions we can prove that the 
local attractors found in 1171 are, indeed, global attractors. 
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we have that for 6/5 < 7, N — > 00 or V — ► 1. Second, a "local analysis" shows 
that the conjecture is true locally. And third, in our numerical analysis no other 
behaviour has been seen. In Appendix IT1 some of the numerical plots are presented; 
in particular, we see that M = 1/N — > for various values of 7. 

3.1. Equilibrium points. 

3.1.1. B{T): equilibrium points of Bianchi type I. 

(1) 1(1): q = i(3 7 - 2), Q = 1, £+ = S x = S x = N x = N = V = 
Eigenvalues: ±(37~2)[x3], -±(2- 7 )[x4]. 
The remaining equilibrium points of type I are all in dT(I). 

3.1.2. T(II): equilibrium points of Bianchi type II. All of the equilibrium points 
in the set T(II) are unstable. They are all given in j^U] (in 'F-gauge'). 

3.1.3. T(VIIq): equilibrium points of Bianchi type VIIq. Using the monotonic func- 
tions Z 3 and Z 2 , it is possible to show that there are no equilibrium points in 
T(VIIo), apart from the following: 

(1) Cl{VII ): q = 2, E + = -1, N = k = constant, S x = Si = N x = ft = 0, 
with the tilt velocities given by: 

(a) < 7 < 2, V! = v = 0. 

(b) < 7 < 2, vi = 1, v = 0. 

(c) < 7 < 2, vi = 0, |v| = 1. 

(d) 7=1, vi =0, 0< |v| < 1. 

(2) T k {VIh): q = 0, N = k = constant, = 1, E + = S x = S x = N x = V = 



Eigenvalues: 0, -2[x2], -1 ± sf\ - 12fc 2 [x2]. 
This is an attracting set (for 7 = 2/3). 

All of the equilibrium points in T(VIIq), including the ones in the sets T(II), B(I), 
dT(I), can be shown to be unstable into the future for 2/3 < 7 < 2. 

3.2. Late time analysis. For the type VIIo model it is convenient to solve con- 
straint (|15(l to obtain an expression for Si: 





In the following we will introduce the variable 



(22) 





where B is a bounded (complex) function. 
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Due to the oscillatory behaviour of the system, we introduce the following vari- 
ables 



(24) 
(25) 
(26) 

where 

(27) </>' 

(28) 9' 



iN x = e 2M/, X, 
iN x = e~ 2 **Y, 
v = e- l9 v 2 , 



M 
M 



1 



Vl 



-Mlm 



(l + wi)Xe 2i ( e+ ^ + (l-t-i)Ye 



2i(0-i/>) 



The angular variables tp and 9 are introduced to take care of the rapid oscillation as 
M — » 0. We note that the variables ip and 9 are not in synchronization since v\ < 1. 
Hence, in general, we expect two different oscillations with different frequencies. 
The equation of motions for these variables are given in Appendix El In this case 
we can explicitly see how the oscillatory terms enter into the equations of motion. 
Moreover, we note that both of these rapid oscillations are observable; e.g., by 
considering the scalars: 



Si = 

s 2 = 



(E x +«N X )(E* X + *N* x ) = e 4 ^XY*, 
(£ 2 x +N 2 x )(v*) 4 = e 4 * e XY U2 4 . 



3.3. Reduced System. The idea is that as M — > at late times, the system of 
equations effectively reduces to a much simpler system of equations. In Appendix 
lAl the idea behind this reduction is explained in more detail; by introducing new 
variables this reduction of the full system is manifest as M — > 0. This method was 
also used in the previous analyses of the type VIIo model |21l 1221 fT5| with success. 
Moreover, in Appendix iBl a linear analysis is given to put bounds on solutions of 
the reduced system with respect to the full system. 

Therefore we assume that M — > 0, which implies that Ei — > 0. Consequently, 
from above we have that all oscillatory terms effectively cancel. Defining 

(29) o~\ = |N x | 2 + |S x | 2 = i(|X| 2 + |Y| 2 ), 

(30) a 2 = 2Im(E x N* x ) =-i(|X| 2 -|Y| 2 ), 

the system given by eas. (|83l89|l effectively reduces to the following system: 

(31) S' + 

(32) a[ 

(33) M' 

(34) n' 

(35) v[ 

(36) v' 2 



= (g_2)E + -a 1 + ^-(-2t; 2 + W 2 ) 

= 2(Q + S+-l)ai 

= -(Q + 2S + )M 

= §^{ 2 Q- ( 3 7 - 2) + [2Q( 7 - 1) - (2 - 7 ) - 7 S] V 2 } 

= (T + 2S+) Wl 

= (T-Z + )v 2 
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where 
(37) 
(38) 



= 2Ej 



l (3 7 -2) + (2- 7 )y 2 

CTl H ; : ^ il, 

2 1 + (7-l)l/ 2 



V 2 S = (-2t> 2 + vl) S+. 
These variables are subject to the constraint 
(39) 1 = E 2 , + a x + Q. 

Furthermore, 02 is determined from 



(To = 



which gives the bound 
(40) 



01 > 



7»M 
G + ■ 



Using the constraint (|39(l we can solve for <j\ or f2. 

3.3.1. Monotonic functions for the reduced system. By the same trick as in Appen- 
dix^ monotonic functions of the reduced system will at sufficiently late times also 
be monotonic for the full system. It is therefore useful to list some of the monotonic 
functions for the reduced system. 



20 m= 1(4-37) 



1 — mS 



+ 



2 m(l — m)(2 — m)|v| 2 f2 



Ru 



which is monotonically increasing for 2/3<7<4/3ifv^0, and for 2/3 < 7 < 2 
if v = 0. 



R 2 = 



(1 _F2)|(2- 7 )' 

3(3 7 -4)i? 2 , 



which is monotonically decreasing or increasing for 7 < 4/3 and 7 > 4/3, respec- 
tively. 



(41) 



R3 

R', 



2S + + ii-M[(l-^)(l-fi)+7U 2 ] 



R3, 



which is monotonically increasing for 7 < 4/3. 



3.3.2. Equilibrium points for the reduced system. 

. P\: E+ = <ri = Af = «i = «a = 0, = 1, Q = |(37 - 2), 2/3 < 7 < 2. 
Eigenvalues: 

™(37-2), -(4-3 7 )[x3], -|(2-7). 
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P". V - 1 Cir V A\ „ - 3(2- 7 )[(3 7 -4)-^(5 7 -4)] 

■ ^2- h + - - 5 W - 4J, <7i — 4(1-^) : 

Q = |(37-2), < i>i < 4/3 < 7 < 2. 
Eigenvalues: 



m = , 2 = o, n = 



-3(4-37), -|(2-7), 



4[(3 7 -4)+^(5 7 -4)] 
V 1"(7-1)«? 



. P 3 : S + - M = 0, ax = ffjff^ 
(2 + 3wi)«i < wf +u| < 1. 
Eigenvalues: 



= 



_ 3+«£-j-v| 



3(1-^1)' = X > 7 = s, < 



0fx2l 



' 2 {1 - V f +v 2 ){3 _ V 2 ) J 



. P 4 : £+ = M = 0,ax = ^fy, «2 = V 1 - « = 3(1=^), Q = 1, < i>i < 1/3, 

2/3<7<2. 

Eigenvalues: 



0, 



2(3 7 - 4) 
2-7 



■1, 



l±iJ3 + 12v( 



4. Late time behaviour 



Theorem 4.1. For all tilted Bianchi models of type VIIo with a perfect fluid stiff er 
than radiation (4/3 < 7 < 2) and where the fluid has non-zero vorticity (\ =/= 0), 
the fluid will asymptotically approach a state of extreme tilt; i.e., Iinx,-—^ V = 1. 

Proof. From the monotonic function Z3 we have that V — > 1 or N — ► 00. First, 
we assume that N — > 00. At sufficiently late times we can then use the reduced 
system. For v\ ^ 0, we use the function R2 which is monotonically increasing at 
sufficiently late times. This shows that V — ► 1. We assume, therefore, that vx = 0. 
For vx — there is a monotonic function given by 



Ra 



ax v 



(l_ya)(2- 7 )> 



1 - y 2 

(37-4) ( 1 + ^^fi 



i?4, 



which again implies V — ► 1. The theorem now follows. 



□ 



As summarized in Tabled the tilt becomes extreme only in the case | < 7 < 2 
for fluid with vorticity. In the case of zero vorticity, a very different tilt behaviour 
is observed; namely, the tilt tends to a non-extreme limit, as described by the 
equilibrium point Pi |19l 125) . 

Theorem 4.2. Consider a perfect fluid Bianchi type VIIo model with 2/3 < 7 < 
4/3 for which N — > 00. IfV < 1, then: 

(1) lim T ^ 00 (£+, V, ax) = (0, 0, 0) for 2/3 < 7 < 4/3, 



(2) lim r _ 



for 7 = 4/3. 



Proof. Use R3. 



a 
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4.1. Decay rates for the case | < 7 < |. By assuming an ansatz with coeffi- 
cients and exponents to be determined, we obtain the decay rates when 3 < 7 < § 



as follows: 






(42) 


£+p 




(43) 


£71 P 


d £7ie-( 4 - 3 ^, 


(44) 


Vi P 


d ie -( 4 -37)T ) 


(45) 






(46) 


M p 





The angular variables are given asymptotically by: 



(47) -0 ps -0 H -e' 

(3 7 -2)M 



2/3 < 7 < 10/9 



(48) #«<j0+^r, 7 =10/9 

+ (iStM el(97 - 10)T ' 10 / 9 < 1 < 4/3- 

The meaning of the bifurcation value 7 = 10/9 is more clear if we calculate the 
fluid vorticity; in particular, to leading order we have: 

(49) W a W a oc e (97 - 10)T , W°W oc e 3 ^" 6 ^. 

We note that the bifurcation values for the vorticity and tilt at 7 = 10/9 and 
7 = 4/3, respectively, coincide with with the values found in |26j - 
The Hubble-normalised Weyl scalars Wi and W 2 , defined by 

(50) m = afccrf w? = bcd 

[ ' 1 48tf 4 ' 2 48tf 4 ' 



evolve as 



12 



(51) => (Wi,W 2 ) « -^- v /^^2e 6(7 - 1)T (cos20 2 ,sin20 2 ) , 

where < Vi < — , M is the constant from the variable M, and 02 is defined by 
Si = \Si |e 2 "^ 2 . The angular variable -02 is related asymptotically to the variable 
via 02 ~ 20 + 0o, where 0o is a constant. 

We note that the first term in 142f) is dominant, and the second term is included 
to display the constant E+. The Weyl scalars diverge for 1 < 7 < |. In the case 
7=1, the magnitude of the Weyl scalars are asymptotically constant. 

The decay rates are verified by numerical simulations. 



4.2. Decay rates for the radiation case: 7 = |. The quantitative late-time 
dynamics for the case of radiation (7 = |) is of particular interest. The reduced 
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system in this case is 

(52) £' + = -S + + E 3 + - ai + \ 2 (1 - Si - <n)(-2w? + u|) 



(53) a[ = 2(1 + 
6(1 

"i 



( 54 ) «i = o 2 :.2 s +^i 

O V-i 



/r^ / 3(1 + «? — 1^) _ 

(55) 4 = i- ± ^S+i; 2 



with the bound 



4(1-2=-^), 



(56) a x > ^ ■ 2 + , 2 > i 



1 ^ "2 



In view of the limit P3, we assume the ansatz 



(57) 


£+ = 


--E + e 


-fer 




(58) 


Cl = 


= a\ + 


(Tie" 


-feT 


(59) 


Vi - 


= V\ + 




-fer 


(60) 


V2 = 


= #2 + 


«2e" 





where fc, the hatted and the barred variables are constants. The ansatz is substi- 
tuted into the evolution equations and terms with equal power are matched. We 
determine the following constants: 

2{vl - 2vf) 



(61) ax 
(62) 
(63) 

(64) 



3(1 -v 2 x+v 2 2 ) 
2S+0-1 




That there are two solutions for k means that the ansatz should have been extended 
to include two modes. The remaining arbitrary constants are vx , V2 and two £+(±) 's. 
The bound lf56*|l restricts vx and £2 as follows: 

(66) 2vx(l+vx) <v% < l-v\ , 0<ui<i. 

We note that when the k± are complex, the late-time approach is oscillatory. 
The angular variables are given asymptotically by: 

(67) 
(68) 

M 
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The Weyl scalars diverge as 



. . /Cf,2 _ 7 -,2\2 _ 4f.2 

(69) (Wi , W 2 ) « - Vla 1 e 2 ^ (cos 2</>2 , sin 2^ 2 ) , 



where M is the constant from M ~ Me T , and ^2 is again defined by Si = \Si \e 2t ^ 2 . 

4.3. Decay rates for the case | < 7 < 2. 

(70) M w Afe~ T 

„ 2(3-y— 4) 

(71) 1 - V 2 « C (1 _ V 2 )e -^^ r 



and E + , tri, «i and V2 have the form 
(72) 



y w j/ + C 1 e~"^r ±T + C 2 e-i T cos 3 + 12vf t) + C z eT* r sin \^J3 + \2v\ rj , 
with 

(73) £ + =0, = 3 1 ( ~_ 3 | ) , v 2 = ^JT^. 

The actual expressions for the constants are not important. The free constants are 
M, C(i_y2 ); C s+ 2 and C s+3 . 

The angular variables are given asymptotically by: 

(74) 
(75) 

M 

The Weyl scalars diverge as 



12 / 1 — % 2 

(76) ( Wl ; w 2 ) « _ . / 1 e 2r (cos 2 ^ 2 , sin 2V-2 



M 2 V 9(1 - v\) 



where < v\ < |. 



5. Discussion 

We have analysed the asymptotic dynamical behaviour of the tilted Bianchi 
type VIIo model at late times. The analytical results are supported by numerical 
calculations (which are discussed in Appendix C). A summary of the late time 
asymptotic behaviour of the tilted Bianchi type VIIo model is given in Tabled The 
attractors all refer to the reduced system, which is valid in the limit M = 1/N — ► 0. 
A striking feature of the type VIIo model is the behaviour of the Weyl tensor. For 
1 < 7, the expansion-normalised Weyl scalars Wi and W2 are unbounded into the 
future. In the terminology of |27| . the model is said to be extremely Weyl dominant 
at late times. Note that the (non-expansion-normalised) Weyl invariants actually 
decay at late times; e.g., for 7 > 4/3 

(C abcd C abcd , C abcd *C abcd )=A%H\W u W 2 )~e- 6T (cos2V>2, sin2^ 2 ). 
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Invariant 










subspace 


Matter 


Attractor 


Comments 




2/3 < v < 4/3 




Pi 


non-tilted 




7 = 4/3 


P 3 


vortic 




4/3 < 7 < 2 




P 4 


V — > 1, vortic 


F(VII ) 


2/3 < 7 < 4/3 


Pi 


non-tilted 




7-4/3 


P 3 


Vi=Q 


vortic 




4/3 < 7 < 2 


Pi 


Ui=0 


V — ► 1, vortic 


T^VIIo) 


2/3 < 7 < 4/3 


Pi 


non-tilted 




4/3 < 7 < 2 




P 2 


tilted 


B{VII ) 


2/3 < 7 < 4/3 


Pi 






4/3 < 7 < 2 


P2 


vi=Q 





Table 1 . The late-time behaviour of the VIIo Bianchi model with 
a tilted 7-law perfect fluid (see the text for details and references). 
The comments refer to the late-time asymptotics, and for all cases 
N — > 00. The case < 7 < 2/3 is covered by the no-hair theorem 
(the non-tilted version is given in [31], the tilted version in 
The value 7 = 2/3 is a bifurcation value for which N — ► constant. 



The general tilted type VIII model, which shares many of the features of the 
Bianchi type VIIo model with a tilted fluid studied here, is currently under in- 
vestigation. In order to complete the analysis of the late-time behaviour of tilted 
Bianchi models, it then remains to study the Bianchi type VI^ models 5 . 

We have focused on the late-time behaviour of the tilted Bianchi type VIIo model. 
However, since the Bianchi type VIIo model has type II as part of its boundary, it 
is plausible that the tilted type VIIo model is chaotic [3131 at early times close to 
the initial singularity. 

In this paper we have utilized a formalism adapted to the timelike geodesies 
orthogonal to the hypersurfaces defined by the type VIIo group action. We note 
that if we consider the integral 

/•oo 1 roc g _ 5 

(77) AT= / —y/l- V 2 dr ~ / e^T r dr, 

which corresponds to the proper time of an observer from t = tq to r = 00 following 
the fluid congurence, there is a change of behaviour at 7 = 8/5. For models with 
7 < 8/5, we find that this integral diverges. However, for fluids with 7 > 8/5, 
we find that this integral is finite. Since asymptotically at late times a fluid with 
7 > 8/5 will be extremely tilted, this means that the fluid will reach future null 
infinity in finite proper time. So in spite of the fact that these spacetimes are future 
geodesically complete the world-lines defined by the fluid congruences seem to 
approach the boundary so quickly that they reach null infinity within finite proper 
time as measured by the fluid. A similiar phenomenon occurs for the LRS type V 
model with 7 > 4/3 The physical interpretation of these models is extremely 



We should stress that some results regarding the stability of certain solutions are already 
known ISBIIT5II25| . 
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complicated. To fully understand this behaviour we need to study the dynamics 
using a formulation adapted to the fluid (i.e., a fluid- comoving viewpoint). We shall 
return to this in future work. 
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Appendix A. New variables 
Let us introduce the following variables 



(78) 
(79) 
(80) 
where 

(81) V' 

(82) 6' 



•»N X = 
*N X = 
v = 



= e 2 ^ 
e 



a/3 | Nv! 



-Im 



(l+^)Xe 2, < 9+ * + (l- Wl )Ye 



From the remaining freedom we have in choosing the variables and the gauge func- 
tion we can, for example, choose the initial values for X and Y to both be real. In 
any case, objects like |X| 2 and |Y| 2 are gauge-independent and are consequently 
independent of any such choice. We also define M = 1/N. 



K 

(83) 

X' 

(84) 

Y' 

(85) 
(86) M' 
fl' 

(87) 

v[ 

(88) 
(89) v' 2 



(Q 2)£ + + 3| Sl | 2 - \ (|X| 2 + |Y| 2 ) + Jg- (-2vj + vt) 

+ (£+ + l)Re (X*Ye- 4 ^) , 
(Q + £+ - 1)X + V3S 2 e- 2 ^ 



+Rc (X*Ye- 4 ^) X - (1 + £ + )Ye _4i * + 



2G+ 



(Q + £+ - 1)Y + V3S 2 e 2 ^ 
+Re (X* Ye- 4s *) Y - (1 + E+)Xe 4 ^ + 

-M [Q + 2S+ + Re (X*Ye" 4 ^) + V3MIm (Y'Xe 4 **)] 
-^-{2Q - (37 - 2) + [2Q( 7 - 1) - (2 - 7) - 7< S] V 2 } 

+2Rc (X*Ye-^) ft 
(T + 2S+) wi - 2V3Re (£ie* e ) v 2 
. V3 



Re (xe 2i < e+ ^ - Ye 2l < e -«) v\ 

(1 + Vl )Xe 2i ( d+ ^ + (1 - Vl )Ye^ e -^ 



T-£+-^Re 



V2, 
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where 



(90) Q - + IRf) + 1 0X|- + |Y| 2 ) + 1 (37 1 - + 2) (7 + _ <2 1 - v 7f ! 
These variables are subject to the constraints 

(91) 1 = £2 + | Sl |2 + l(| X |2 + | Y |2) +fi 

(92) = I(|X| 2 -|Y| 2 )-2gl. 

We can also split V 2 S and T into oscillatory, and non-oscillatory parts: 
V 2 S = H+(-2vl + vl) 

(3 7 - 4)(1 - V 2 ) + (2 - 7 )£+ (-2v\ + vl) 



2 



(93) +^Re 



T : l-( 7 -l)V 2 

(94) + 2(1^(7-1 7 )V 2 ) RC [(^^^ + Ye2l(e ^') «2 + 4S 1 e i V*, 2 ' 

A.l. Assuming M — > 0. Let us assume that M — > in the following. To inves- 
tigate the asymptotic dynamics we choose a new set of variables. If x is a generic 
variable, we define x as: 

(95) , = _£L__M 5 , 

where / and g are bounded functions of the state space variables. If the equation 
of motion for x is 

(96) x' = Fx + G, F,G bounded, 
we obtain 

-/ (F -Mf)x + G-Mg' „,„ 

(97) x 1 = * 1. + MB, 

where B is a bounded function into the future. The trick now is to choose the 
functions / and g in order to get rid of the oscillatory terms and write the system 
in the following form: 

(98) x' = Fx + G + MB, 

where the function B is bounded to the future. For the variable M, we choose the 
following form: 

(99) M' = FM + M 2 B. 
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Using the transformation ea . I|95|) . the following functions / and g will do the 
trick for the various variables: 

£+ : / = g = ^Im (X*Y e - 4 "^ 



X : 

Y : 

M : 
Q : 

Vi : 



/ 
9 : 
f 
9 ■■ 
f 



1 



4V3 

i 

4\/3 
1 



4a/3 

Im (X*Ye- 4 ^) , 
1 + S + )Ye- 4 ^ 
Im fX*Ye- 4 ^) , 



4G+(1 + Ui) 



2 

2 e _2i(0+VO 



4^3 
1 

I7s 



( 1 + ^ X ^'-4G + (l- Vl) 
Im(X*Ye- 4 ^), <? = 0, 



2V3 



Im fX*Ye" 4 ^) 



4G 4 



-Im 



,-2i(e-v) 



_2^|_ e 2 l (e+v) 



1 + Ul 



1 — Ul 



5 = 0, 



«2 



/ = 

9 = 
f = 



(2-7) 
4G_ 
1 

Im 

4 

(2-7) 
4G_ 



111! 



_2^2_ e 2i(9+V) _ Y "2 e 2i(6-Tp) 



1 + Vi 



l- Vl 



1 + TJl 1 — Ul 



Im 



_2^2_ e 2i(0+VO _ Yi, 2 e 2i(0-Tp) 



Im 

4 



1 + v\ 1 — V\ 

5 = 0. 

We note that all of the above functions / and g are bounded; in particular, 



Xe 2i(<?+V>) „ Ye 2i(e 



< 



< 2. 



l±«i 

Let /i* be the transformation /i l : a; 1 1— > i 1 given above (supplemented with 
(9,ip) — (9,ip)). We note that if there exists an e > such that < 1 — e, then 
the Jacobian J = 



fats') 



has determinant 

k 



(100) 



det(J) = 1 + ^2 Mn bn, b n bounded. 



71= 1 



Hence, by the inverse function theorem, if M — > there will exist a to such that 
h % : x 1 1— > i l has a continuous inverse for all t > to. This means that close to 
the equilibrium points of the reduced system, the functions h l and their inverse 
are well defined. We can consequently rewrite the system as ea. (|9"5j) close to the 
reduced system. As M — ► 0, we can then use centre manifold theory to show that 
the system effectively reduces to the reduced system of section 



Appendix B. Bounds on oscillatory terms 



Consider an equilibrium point, Xo, for the reduced system. If we linearise the 
full system about X (treating 9 and ijj as sources and not part of the system) we 
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get a matrix equation of the form 

(101) X' = (L + S(0,V))X + P(0,V). 

Here, the matrix S(0,ip) and the vector P(0,ip) contains the oscillating terms that 
need to be estimated. 

This linearised equation has the formal solution 

(102) X(i 2 ) = e^'^Xih) + e I(tl ' t2) [*' e" 1 ^ 1 '*^^, if>)dt, 

Jti 



where 



I(h,t 2 )= f 2 (L + S(0,V))di. 
Ju 



In order to estimate the contribution from the oscillating terms we need the 
following result: 

Lemma B.l. Consider the integral 

(103) i{h,t 2 ) = f Vv^dt, p,p>0, p>\. 

Then there exists a constant D such that for every t\ and t 2 satisfying 

h < t 2 , 

we have 

(104) \I{h,t 2 )\ < D (V A - p)tl + e (A -" )t2 ) . 

Proof. We introduce the variable z — e pt so that the integral can be written 

1 f Z2 » -a 
I = - z—e l0z dz. 

P J Z1 

Considering z as a complex variable, we notice that the integrand only has (at 
most) a singularity at z = 0. Hence, the integrand will be analytic in the entire 
half-plane Rc(z) > 0. By choosing a suitable closed path we can rewrite I as: 

(105) I = -7 Cl - I C2 , Ci = Zi + iy, 0<y < oo. 
These two integrals can be estimated: 



(106) p\I c 



f> OO /'OO 

/ ( Zl +iy)^ E e l()zi - f3y dy < / (zf + y 2 )^ e'^dy. 
Jo Jo 



Since A < p the factor [z\ + y 2 ) 2 " is monotonically decreasing for y > 0; hence, 

a-p r°° a i a-p 

(107) p\I Cl \ <z x " ^ e-^dy=-z^ . 

Thus we have 

(108) |/| < |/ Ci | + |7 C2 |<-L^+^^ = _L( e (A-^ 1+e (A-^ 2 ^ 
which proves the Lemma. □ 
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Consider first the linearised M equation with respect to the future stable points 
of the reduced equation. We assume that M — > 0, and by a similar complex 
integration we get 

{109)M(t 2 ) = e-^- t ^M(t 1 )+O(M 2 ), p =\^ 3nf ~ 2 ^ f < 7 < 3 

[1, 3^ <2 ' 

Hence, asymptotically, M oc e~ pt . 

The matrix S(9,ip) contains oscillating terms like e M , where a is a linear com- 
bination of 9 and if). Using the equation for 9 and ip, and the leading order approx- 
imation for M, we obtain 

(110) I(ti,t 2 )= / 2 (L + S(M))dt = L(t 2 -«i) +0(e-" tl + e-" t2 ). 

Let us expand X (and P) in terms of a sum of eigenvectors X\ of L; i.e. 

X = ^Xa, LXa = AXa- 

A 

The solution can now be estimated: 

(111) X A (i a ) = e A (* 2 -* l) X A (ii) + o( e -p*i+ A (*2-*i) + e ~pt2y 

This result shows that the deviation from the solution of the reduced system is of 
the order of M. Hence, since M — *■ 0, the linearised solution will asymptotically 
approach the solution of the reduced system. 

Appendix C. Numerical analysis 

Numerical calculations can be used to both confirm and clarify the analytical 
calculations found in the bulk of the paper. We have chosen to use the new vari- 
ables found in Appendix A, in which the oscillatory part of the system of differential 
equations can essentially be isolated from the non-oscillatory part through the in- 
troduction of the variables tp and 9. Equations H81(l - (|89[1 (minus equation 187(1 ~) 
were numerically integrated using the 'ODE23t' ODE solver in Matlab with an 
Absolute Tolerance of 10~ 6 and a Relative Error of 10~ 6 . The constraint equation 
11921 was used to determine if and when the numerical calculations broke down. 

Many different sets of initial conditions were tested to confirm the results of the 
analytical calculations. We observed that M — * in every numerical run for values 
of 2/3 < 7 < 2. We also observed that the final asymptotic state obtained in the 
numerical calculations confirmed what was found in the analytical calculations. In 
Figures 131 71 we show some of the numerical runs for different representative values 
of 7. The value 7 = 19/15 is chosen since it represents a typical value between the 
'special' value of 7 = 6/5 and the bifurcation value 7 = 4/3. 
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Tilted Bianchi VII (y=1) 



Tilted Bianchi VII (7=1 ) 



10 12 




10 12 



FIGURE 3. Plots of M, £+, v\ and |X| 2 for 7 = 1, illustrating the 
rate at which M (and all other variables) tend to zero. 
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2 

°D 
1 

0.8 
0.6 
0.4 
0.2 

0, 



Tilted Bianchi VII (y=19/15) 



2 4 6 8 10 
X 

Tilted Bianchi VII (y=19/15) 



Tilted Bianchi VII (y=19/15) 



0.8 
0.6 
0.4 

W + 0.2 


-0.2 
-0.4. 




10 12 



Tilted Bianchi VII (7=1 9/1 5) 




FIGURE 4. Plots of M, S + , v\ and |X| 2 for 7 = 19/15, illustrating 
the rate at which M tends to zero. Compare the slow approach of 
£_!_,«! and |X| 2 towards zero with the corresponding plots with 
7 = 1. 
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Tilted Bianchi VI l Q (y=4/3) 



Tilted Bianchi VI l Q (y=4/3) 




6 8 10 
X 

Tilted Bianchi VII Q (y=4/3) 
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0.8 
0.6 
0.4 
0.2 





1 

0.8 
0.6 
"0.4 
0.2 




Tilted Bianchi Vll (y=4/3) 




Figure 5. Plots of M, £ + , vi and |X| 2 for 7 = 4/3, illustrating 
the rate at which M tends to zero. Note how the variables v\ and 
|X| 2 do not approach a unique value. 



Tilted Bianchi VI l Q (y=5/3) 



Tilted Bianchi VI l Q (y=5/3) 




4 6 8 10 
X 

Tilted Bianchi VI l Q (y=5/3) 




Tilted Bianchi Vll (y=5/3) 




FIGURE 6. Plots of M, S + , vi and |X| 2 for 7 = 5/3, illustrating 
the rate at which M tends to zero. Note how the variables v\ and 
|X| 2 do not approach a unique value. Also notice the oscillatory 
behavior and how it dampens. 
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Figure 7. Phase portraits of v% versus «2 for different values of 
7 illustrating the fundamental differences in the asymptotic states 
in each case. Note how v\ + v\ ^ for 7=1 and 7 = 19/15, 
v\ + v\ -> y 2 ^ for 7 = 4/3 and v\ + v\ -> 1 for 7 = 5/3. 
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